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Dedicated to the memory of Professor Leon Ehrenpreis, 
a great mathematician and human being. 

Abstract. The article describes interesting nonlinear differential 
identities satisfied by standard exponential and trigonometric func- 
tions, which appeared as byproducts of medical imaging research. 
They look like some kind of non-commutative binomial formulas. 
A brief description of the origin of these identities is provided, as 
well as their direct algebraic derivation. Relations with separate 
analyticity theorems in several complex variables and some open 
problems are also mentioned. 



1. Introduction 

Working quite some time ago on some problems of medical imag- 
ing [I9l|20j, the authors were surprised when their work produced 
a series of seemingly new (at least to them) identities for functions 
about which "everything is known," such as standard exponentials and 
trigonometric and hyperbolic sines and cosines. It was a little bit de- 
pressing when our brave attempt to prove these identities by elementary 
means, without using the mathematical techniques common in tomog- 
raphy (c'mon, what can be difficult about sinx, after all?), stalled for 
a while. We have succeeded eventually, albeit as we hope to persuade 
the reader, there are still things that are not clear about this whole 
business. Along the way, while trying to understand the identities, an 
interesting theorem of complex analysis in several variables was discov- 
ered in the joint work [Tj of one of the authors and his student with L. 
Ehrenpreis. We dedicate this article to the memory of Professor Leon 
Ehrenpreis, who passed away in August 2010 (see pniB3]). 

In the next Section [2] we describe briefly the medical imaging origin 
of the results we will be discussing. Section [3] contains the formulation 
of the identities, which are then proved in the next three Sections. 
The paper ends with the remarks section, where in particular some 
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open questions and relations to several complex variables theory are 
described. 

2. A BRIEF STORY ABOUT MEDICAL IMAGING, RADON 
TRANSFORM, AND ALL THAT 

In tomographic medical imaging one sends through the patient's 
body some kind of radiation/ waves (X-rays, ultrasound waves, elec- 
tromagnetic waves, etc.), measures the transmitted and/or reflected 
signals, and tries to recover from this data the internal distribution of 
some important for medical diagnostics internal parameter of the body 
(e.g., density, electrical conductivity, blood oxygenation, etc.). 

The reader must have heard of at least some tomographic proce- 
dures (the best-known ones are the X-ray CT scan and MRl). How- 
ever, not everyone might realize the heavy involvement of sophisti- 
cated, beautiful, and challenging mathematics in creating tomographic 
images. In spite of its 40 years history, computed tomography still 
produces a steady stream of exciting mathematical problems (see, e.g., 
[ai5l[71|9l[IOl[IlMl[2ll|2M5l[28l[29l|3^ for description and math- 

ematics of established and newly developing imaging modalities). 

In the oldest and best established modality. X-ray CT, one deals 
with the problem of recovering a function f{x) that is, roughly, the 
tissue density at the internal point x. Since in many cases different tis- 
sues have different densities, this function carries medically important 
information about the interior of the patient's body. If f{x) is found, 
its density plot provides a picture, tomogram that shows to a doctor 
how an internal slice of the patient's body "looks like." Mathemati- 
cally speaking, the data obtained by the X-ray CT scanner provides 
the practitioner with line integrals of the sought for function / (see, 
e.g., [7l[ini[IZl[22l[21] for detailed explanation). In other words, the 
data one measures is the so-called 2D Radon transform that takes 
function f{x) on the plan^ and integrates it over all lines L: 



We will parametrize lines by their normal coordinates {u, s), where 
w is a unit vector normal to L and s is the (signed) distance from the 
origin to L. Then the normal equation of the line L has the form 



We will not be precise about conditions on functions we deal with here. It is 
safe to assume that they are "very nice" : differentiable as many times as you want 
and vanishing at large distances. 




L 



X ■ OJ = s 
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where x ■ uj is the standard inner product. Thus, the Radon transform 
Rf of a function / can be written as follows: 



(1) Rf{uj,s) = g{oo,s):= / f{x)dl = f{soj + toj^)dt 



Here a;-*- = {u2, —oJi) is the rotated 90° vector u. 

The Radon transform, besides being extremely useful in tomography, 
arises in a large number of other areas, e.g. in PDEs, group representa- 
tion theory, etc., and thus has been the topic of thorough mathematical 

study [Sl[I2HISl[2a[23[2Sll2S]. 

We are interested here in one feature particular for the Radon trans- 
form: its range as of a linear operator between natural function spaces 
is rather small, i.e. of infinite co-dimension. In other words, there 
are infinitely many conditions a function g{(jj, s) must satisfy to be the 
Radon transform of a function /. Let us discuss them briefly. 

One such condition is evenness: 

(2) g{—<^, —s) = g{u), s) for all uj G SS, s G M, 

where SS is the unit circle. This reflects the fact that equations x-uj = s 
and X ■ (— = —s describe the same line. 

More interesting are the so-called moment conditions. To formulate 
them, we assume that / has a compact support (in fact, decay faster 
than any power of |x| at infinity suffices). Then for any integer k > 
one can define the k-th moment of the data g{u,s): 



(3) Gkico):= / s'g{LO,s)ds 



which is a function on the circle SS. 
Now the moment conditions say that 



For any integer A; > 0, the moment Gk{uj) 
(4) can be extended from SS to the whole plane 

as a homogeneous polynomial of degree k. 

Although it might not be immediately obvious how one discovers the 
moment conditions, as soon as they are formulated, it is straightforward 
to check their necessity. Indeed, substituting into ([3]) the definition of 
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the Radon transform g = Rf, one gets 

OO OO 

Gk{u}) := / s''g{uj,s)ds = s'' f{x)dlds= I {x ■ lo)'' f{x)dx. 



The last expression is clearly a homogeneous polynomial of degree k of 
u (since x ■ u is a homogeneous linear function of u) and thus we have 
established these so-called "range conditions" . 

Proving that there are no other range conditions besides evenness 
and moment ones, is a different story, which takes much more work 
(see, e.g. [5l[l3l[Tl[2a[25j ) . 

Why would one be interested in these conditions? They are the 
mandatory relations that ideal data collected from a tomographic de- 
vice must satisfy. Well, the measured data are never ideal, and so 
they deviate from these conditions. Thus, knowing the range condi- 
tions might be helpful in detecting and correcting some measurement 
errors. They are of help in other circumstances as well, for instance in 
completing some missing data (e.g., see [22l[21l[2SlEn] and references 
therein for details). And they surely play an important role in most of 
analysis of the Radon transform as an operator. 

Quite a few years ago, the authors worked [T9|l20] on finding the range 
conditions for some special (weighted) Radon-type transform arising 
in another popular medical imaging method, so-called SPECT (Single 
Photon Emission Computed Tomography) [9l[T6 |[T8l[2n[2l] . This is the 
so-called exponential Radon transform, which is defined as follows: 

OO 

(5) RJ{cu,s)=g{cu,s):= j f{sco + tco^)e'''dt. 

— OO 

Here /i > is the attenuation coefficient. 

When we found a set of range condition^ for i?^ and then proved 
that we had not missed any, we thought that checking their necessity 
(as long as we already knew them) should be a piece of cake, just 
like in the example of the usual Radon transform above. You plug 
your transform i?^/ into the conditions and should be able to "see" 
immediately that they are satisfied. Well, when we did this, we saw an 
infinite series of identities for sin a; (yes, the usual sine!) that we could 
not recognize and did not know why they should have been true. Here 
they are: 



^We do not provide these conditions here, since they are rather technical and do 
not make what foUows easier. The interested reader can find them in 
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For any odd natural number n, the following identity holds: 
^ f n \ 

E ^ -sinx)o(- -sinx + ^)o... 



(6) to V ^ 



o 



J- — sinx — (A; — l)i)(sinx)" = 0. 



n 



Here i is the imaginary unit, \^ j is the binomial coefficient "n choose 

/c," the expressions in parentheses are considered as differential opera- 
tors, and o denotes their composition. For instance, 

d . \ , , du 



— — smx Ml XI = — — smx u\x). 
dx J dx 

Although we have provided the necessary proofs and moved on to doing 
other things, the true meaning of these identities has kept us puzzling 
for all these years. 

In this paper we present the formulation and elementary algebraic 
derivation of these identities (not only for sinx, but even for some more 
"elementary" functions such as linear and exponential ones). 

We would be happy if someone could shed some more light onto the 
meaning and possible generalizations of these identities. 

3. Algebraic formulation of the identities 

The identities (E]) can be formulated and proven in a rather general 
algebraic setting. Namely, let A be a commutative algebra with cance- 
lation property and with a differentiatioiiH D over a field K. An element 
u & A can be considered an analog of the sine function if it satisfies 
the equation D^u = for some \ & K. For instance, in the case 
when u = sinx, we are talking about the algebra C°°(M) over complex 
numbers of all smooth functions on the variable x, with differentiation 
D = d/dx and X = i. 

In what follows, the reader will not loose anything thinking, instead 
of the general algebraic situation, of smooth functions /(x) with the 
usual derivative D = d/dx. However, the general algebraic notations 
{A, D, A) lead to somewhat easier to read formulas, so we will stuck 
with them most of the time. 

We are interested in solutions of simple differential equations 

(7) Du = \u 
and 

(8) D\ = \\, 



^Le., D{uv) = D{u)v + uD{v). 
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where X G K. Thus, in the example of the function algebra with usual 
differentiation, solutions of equation ([7]) for A = are constants, while 
for non zero values of A they give us exponential functions. Analo- 
gously, solutions of ([8]) for A = are linear functions, while for non 
zero values of A they in particular give us trigonometric and hyperbolic 
sine and cosine. 

We now formulate the main results. 

Theorem 1. Any solution u of the first- order equation Du = Xu sat- 
isfies for any natural n the following identity: 



(9) ^( J: \{D-u)o{D-u + X)o...o{D-u + {k-l)X)u''-^ = Q. 

Theorem 2. Any solution u of the second-order equation D'^u = X^u 
satisfies for any odd natural n the identity 

The expression u^'^ at the very right is just the {n — k)Wi power 
of u. Expressions like {D — u) are considered as operators on the 
algebra A. For instance, (D — u) applied to an element / G A produces 
[D — u)f = Df — uf . The little circles o mean composition of these 
operators from right to left. For instance, 

(D - u) o (B - 1, + A)/ = - "(^))(^ - u(x)f(x) + Xf(x)) 

So, the expression in Q directs us to raise u to the power {n — k), 
then apply to it the operator {D — u + {k — 1)X), then to apply to 
the resulting function the operator {D — u + {k — 2)X) , etc. To avoid 
confusion in the case of = 0, it might be easier to think that the 
composition in ([9]) has k operator factors. In other words, the term for 
/c = does not have any of these factors and thus is just u^. 

An important thing to notice is that, as Calculus teaches us, the 
operations D and u (differentiation and multiplication by u) do not 
commute, i.e. the result of their composition depends on its order. If 
this were not so, the identities ([2]) would become trivial, as we will see 
a little bit later. 

Since solutions of the equations ([7])- ([8]) are very simple functions 
(constants and exponential, linear, trigonometric or hyperbolic func- 
tions), let us look at some examples to see whether the identities ([9]) 
might be trivial, or at least look familiar. 
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Take (J?!) for A = 0, so the solution u is just a constant C. Then ((91) 
looks as follows: 

^ ^ - = 0. 

fc=0 ^ ^ 

Since acting on any constant gives zero result, this identity boils 
down to 

f I ) (-c)''C'"-'= = (-C + C)" = 0. 

fc=0 ^ ^ 



This is indeed easy and proves the particular case of Theorem [T] for 
A = 0. 

Consider now the particular case of ([7]) when A = 1 and u = e^. 
Then (Ml reduces to 



Hmmm... Does not look familiar? Neither it did to us. Well, we will 
not prove this particular identity, since we are going to prove the more 
general Theorems [1] and [21 

Let us look at an example of a solution of the second-order equation 
D'^u = X^u, say take u = sinx and \ = i. Then the identity dHj) for 
n = 3 becomes (using more familiar ^ instead of D) 
(10) 

sin^ X + 3(- sinx) sin^ x + 3(- sinx) o (- sinx + i) sinx 

d d dx ^ dx 

+ (- sinx) o (- sinx + i) o (- sinx + 2i)l = 0. 

dx dx dx 

Obvious, isn't it? Well, do not despair if you do not see it right away, 
we don't either. 

As you will see, the direct algebraic proofs of Theorems [H and El (at 
least the proofs that we know) are different when A = and when 
A 7^ 0. Certainly, it is not hard to derive the identities for A = from 
the ones for A 7^ by taking the limit when A — )■ 0. We, however, 
would like to see direct algebraic proofs of the identities for both cases 
A = and A 7^ 0. Such proofs are provided below. 



4. Proof of Theorem [H 

We have already proven it above for A = 0. So, let us assume that 
A ^ 0. 
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Suppose that u satisfies the first-order equation Du = \u. We can 
apply a gauge transform to simphfy our identities. Namely, we will 
use the easy to check identity 

(11) u-'^i^D -u)u"^ = D -u + mX. 

We can now use ( fTTl) in each factor of (l9l) to rewrite it as follows: 

(12) [j2(^iy{D-u)ou-']^u- = Q. 

n f n \ 

Here we notice that we deal with an operator binomial V ( , \ A^, 

to\k ) 

where A = {D — u) o . Since only one operator is involved, no 
non-commutativity arises, and thus the usual binomial formula works. 
This reduces identity ^ to 

(13) \{D -u)o + 1]"m" = {Do u-^Yu'' = 0. 

Now an immediate calculation using Du = Xu shows that {Dou^^)^u"' = 
holds. 

Remark 3. A careful reader might object to our calculations in this 
section that involve negative powers of u, as well as positive ones, es- 
pecially in the abstract setting of a commutative differential algebra 
A with cancelation. However, this algebraic problem can be overcome 
(see U^). 

5. Proof of Theorem [2] for A = 

We assume now that A = and thus the function u satisfies the 
equation D'^u = 0. In this case the identities that we intend to 
prove simplify to 

(14) Vd-«)V-'= = 

fc=0 ^ ^ 

for any odd natural n. 

Remark 4. This identity fails for n = 2. 



One can understand fll4p also as follows: 
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The expression in the parentheses resembles the familiar binomial ex- 
pression 

fc=0 ^ ^ 

where a = D — u and h = u. If the operators commuted, then the whole 
thing would have been equal to (a + 6)" = [D — u + m)"1 = D"l = 

^ = 0. Done! 

ax"- 

Hold on, not so fast! First of all, since D and u do not commute, 
the usual binomial formula does not apply. Besides, it looks like we 
obtained the identity ^ for all natural n, which is incorrect, as we 
saw in Remark |H So, let us take a closer look. 

If A = 0, then function u is linear. Let us concentrate on the case 
when u{x) = x. 

Let us first rewrite our identity (ITSl) for the case when u = x: 

(16) Y.(l\{D-xfx--' = Q. 

We will now use another popular gauge transformation trick that 
will enable us to rewrite f|T6l) in a slightly different form. It is based on 
the following simple identity: 

(17) iD~x)''{e'^f{x))=e'^D'f{x). 

si 

Thus, commutation with the function e 2 kills the term —x added to 
the derivative and {D — x)^ becomes D^: 

2 2 

e-¥ o {D-xf oe^ = D^. 
This implies the following 

Lemma 5. The identities [TB^] are equivalent to 




So, now the proof of Theorem [2] for A = boils down to proving f|T8|) 
for odd n. Notice that if D and x commuted, then the expression in 
parentheses in flTSl) would become just {D + a;)", and thus, due to the 
easily verified equality 



(19) 



{D + x)e~'^ = 0, 
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we would immediately obtain the validity of ( fT8|) . However, life is not 
so easy and thus the operators do not commute. (Besides, as we already 
know, (fT8|) is guaranteed for odd values of n only). 

Based on (fT8|) . let us introduce the differential operators 



i.— n \ / 



fc=0 

The punch line is in the following statement: 
Lemma 6. The following recurrence relation holds: 

(20) Pn+2iD,x) = Pn+iiD,x)o{D + x) + {n + l)Pn{D,x). 

Proof of the Lemma. We will compute P„+i(D,a;) o [D + x) and 
show that it coincides with P„+2(-D,x) — {n + l)P„(Z),a;). In order to 
do so, we will use the easy to verify commutation relation 

Dox'^ = x"'oD + mx'^-^ 

and its consequence 

(21) x"" o D = D o - mx'"-^ 

Now, take a deep breath, and ... 
(22) 

Pn+i{D, x) o {D + x) = (j2 (^"^^^^ D'' o x'^+i-'^^ o{D + x) 



Let us now use fl2T]) to commute x^^^~^ with D in the first sum to get 



„n—k 



fr)0\ k=0 V / fc=0 

The first two sums can be rewritten as 
^^^^ 

= E ( t ^ ) ° = ^n+2(/^, X). 

We used here the standard property of the Pascal triangle: 
(n + 2)-l\ / (n + 2)-l\_ / n + 2' 



k-1 ' \ k \ k 
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It only remains to handle the remainder 

fc=0 ^ 

This expression is equal to {n + l)Pn{D, x), proving which only requires 
the knowledge of what binomial coefficients are. This finishes the proof 
of the Lemma. □ 

Corollary 7. For any odd natural n, the following factorization holds 
with some operator Q{D,x): 

(25) Pn = Qo{D + x). 

Proof. Let us prove this by induction. When n = 1, we have Pi = D+x, 
and the statement is obvious. Assume that we have proven it for some 
odd n, i.e. Pn = Q o {D + x). Then the previous Lemma implies that 

Pn+2 = {Pn+l + {n+l)Q)o{D + x). □ 

Remark 8. The factorization Pn = Q o [D + x) fails for n = 2. 

Now the proof of Theorem [2] for A = is immediate. Indeed, for 
any odd n, the left hand side in the identity (ITSl) in question becomes 

P„e 2" = Q o (^D + x)e~^ . Applying ( IT9|l . we conclude that this is 
zero. This finishes the proof of the Theorem for A = and u = x. 

The proof provided for u = x can be easily generalized to an arbitrary 
linear function u = ax + h. Alternatively, one can derive the result for 
any linear function from the one for u = x. □ 



6. Proof of Theorem [2] for A 7^ 

Let us now outline the steps of the proof of Theorem [2] in case A 7^ 0. 

Step 1. If M satisfies the first-order equation Du = Xu (and thus 
certainly D^u = X^u as well), the statement follows from Theorem [1] 

Step 2. Suppose Du = —Xu. Then is true for odd n, because one 
can easily check that "symmetric" terms in the sum that correspond 
to k and n — k cancel out. 

Step 3. Let us now suppose that u is any solution of D^u = X^u. 
If our identities were linear with respect to u (which they are not), 
then the cases when Du = iAw discussed in the two previous steps 
would suffice, since any solution could be expanded into a sum of these 
two. Nonlinearity seems to destroy this idea. However, having nothing 
better in mind, let us still try. Thus, u = v + where Dv = —Xv 
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and Dw = \w (notice that the assumption A 7^ is critical for the 
possibihty of such decomposition)@. 

Step 4. Let us use the same commutation trick (fTT!) 

u~™(L> - u)u"' = D-u + mX 

as before, but using commutation with w rather than u. Then the 

left-hand side of iQ can be written as 

(26) 



fe=0 



J2 [ I j (D - V - w) o ... o {D - V - w + {k - l)X){v + w 
[{D -v)o w-^ - l]^{vw-^ + l)"-^' 



in— A; 



n 




E 

.fc=0 





w 



Step 5. Let us introduce the following operator notations: A 
{D — v) o and B = vw~^. Then the last sum becomes 



\n~k 



If the operators A and B commuted, then according to the binomial for- 

1 k TDn—k 



nd B commuted, then accor 

mula this would boil down to (A+B)"' and thus also to ^ ^ ^ ^ A^B"^ 

The interesting thing is that the latter conclusion holds even without 
commutativity: 

Lemma 9. For any two operators A, B the following equality holds: 



E 



{A-l)\B + l 



•.n—k 



k Tjn—k 



A^B 



Indeed, this equality does not require changing the order of the fac- 
tors. Thus, if it holds for commuting operators, then it does for non- 
commuting ones as well. 

How can this help us with the proof of the theorem? It allows us to 
drop the terms =pl in fl26|) to get 



E 

.fc=0 



-l\n—k 



W 



Now one can reverse Step 4 (undoing the commutations with powers 
of w we have done) and rewrite fl26p as 



j2(l)(D-v)o...o{D-v+{k-l)X), 

k=0 ^ ^ 



n—k 



^This is just standard exponential representation of solutions studied in ODEs 
(or, if you will, the representation of sin a; as combination of e^"). 
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which is equal to zero due to the Step 2, since v solves the first-order 
equation Dv = —Xv. This proves Theorem [2] for A 7^ 0. □ 

7. Remarks and generalizations 

• For the readers with a taste for generalizations, we can mention 
that the identities (Q hold in a much more general situation. 
Namely, u can be assumed to be an element of a commutative 
algebra with cancelation and with differentiation D over a field 
A. Then A should be an element of the field A |19j . 

• Generalizations of some of the identities we discussed are avail- 
able: 

Theorem 10. If u satisfies D^u = 0, then the identity 
(27) ( ) - ° ^"^"""^ = 

fc=0 ^ ^ 

holds for all odd integer n > 1 and even m > 0. 

• Several integral geometry and tomography experts devoted their 
time and effort to trying to understand better the meaning of 
these strange range conditions. This is also what one of the 
authors set out to do with V. Aguilar and L. Ehrenpreis. Sur- 
prisingly, as the result, the identities Q) were related [I] to 
the so-called separate analyticity theorems (Hartogs-Bernstein 
theorems) in several complex variables. For instance, the fol- 
lowing amazing theorem is essentially equivalent to these iden- 
tities |ll[26l[27]: 

Theorem 11. Let Q be a disk in and f{x) he a function in 
the exterior ofQ. Suppose that when restricted to any tangent 
line L to dQ, the function fl^, as a function of one real variable 
extends to an entire function on the complexification of L. Then 
f , as a function on \ i7 extends to an entire function on C^. 

Well, this fact also did not look obvious. Analyticity of / in 
a complex neighborhood of \ D follows from the old (and 
not well known) separate analyticity theorem by S. Bernstein 
(see [2j), however this theorem cannot produce statement about 
/ being an entire function. Thus, since proving the above the- 
orem, a couple of things about it kept bothering us for sev- 
eral years. First of all, this is a several complex variables fact, 
while our proof did not look like a SCV argument at all. Is 
there a truly complex analysis proof? Another, related, ques- 
tion is whether such a theorem can be proven for a different 
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convex body instead of a disk fi? A SCV proof was later pro- 
vided in I2SII2Z], although it was rather complicated and was 
not generalizable (at least, easily) to other convex curves. Leon 
Ehrenpreis has worked out some other examples of convex al- 
gebraic curves (unpublished), but general picture remained un- 
clear. Finally, A. Tumanov presented recently [35] a beautiful 
short proof based on attachment of analytic disks (where Tu- 
manov is a great expert), which works for any strictly convex 
body Q with a mild conditions on the smoothness of its bound- 
ary. 

More discussion of the relations of the range conditions with 
complex analysis can be found for instance in |6| IT8] - [20] . 

• It would be interesting to find an algebraic proof of the identities 
that would work simultaneously for A = as well as for A 7^ 0. 

• It is clear that the identities discussed must be related to spe- 
cial function theory and group representations. It would be 
interesting to understand such relations. 

• The formulations of the Theorems [T] and [2] lead to the natural 
question: what can be said for solutions of the equation of 3rd 
order D^u = X'^u and higher? A natural guess would be that 
the same identities hold, but only for an arithmetic sequence of 
numbers n with the difference equal to 3. It is hard to compute 
these expressions by hand even for small values of n, say for 
n = 4. Ms. E. Rodriguez, a former Masters student of P. K., 
has used the Maple symbolic algebra system to check this con- 
jecture. The result is negative, the natural conjecture fails for 
equations of third order [31]. So, what (if anything) happens 
to solutions of higher order differential equations D'^u = X^ul 
Are they deprived of any such identities? We do not know the 
answer to this question. 
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